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Abstract 

In the paper the Schrodinger equation for quasibound resonance state 
with complex energy is considered. The system of inhomogeneous dif- 
ferential equations is obtained for the real and imaginary parts of wave 
function. On the base of known solution of corresponding homogeneous 
equation, the inhomogeneus system is solved with help of iteration proce- 
dure. The single-particle neutron 2p-state in the Woods - Saxon potential 
is analyzed for ^^C nucleus. 

I. INTRODUCTION 

Many papers, published since 1980 and devoted to investigation of the struc- 
ture of nuclei with A = 13, show that for description of the ground states of 
these nuclei, it is necessary to take the 2p - shell into account in a ground 
state [1 - 17]. Recently nuclear structure calculation appeared in [11], which 
described nuclei with A = 4 — 16 within a full (0 + 2)huj shell model space, 
i.e. configuration mixing, is not only restricted to the Ip - shell, but includes 
contributions from the 2s-, Id-, 2p- and l/-shells as well. These wave functions 
provide a good overall description of the relevant energy spectra. Furthermore, 
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from the analysis of the Ml form factor of electron elastic scattering on the 
nucleus ^^C [3 - 11,14], it is followed that the 2p - shell is the most important 
one among the 2s -, Id -, 2p - and If - shells. The similar conclusion regarding 
2p - shell admixtures was obtained for the processes with pion [12,13,15]: pion 
photoproduction, radiative pion capture on ^^C and pion single charge exchage 



If we choose as an average field the oscillator potential, parameters of which 
are obtained in accordance with the binding energy of nucleon in Ip - shell, 
the binding energy of nucleon in 2p - shell becomes positive. Thereby, at the 
value b = 1.663 fm for oscillator parameter and Ei,ind{^p) = —4.947 MeV for 
the neutron binding energy on Ip - shell of ^^C, we find the binding energy in 
2p - shell Ei,ind{2p) = 25.016 MeV in the same oscillator well. (Somewhat less 
value Ei)ind{2pi/2) = 14.65 MeV was obtained in work [18] from the analysis of 
the spectra ^^C, ^^C and ^^C nuclei.) And if the 2p - state with such positive 
energy is bound in the oscillator potential because of the depth of oscillator 
well being infinity, then for the realistic potential of Woods - Saxon, the 2p - 
state represents itself as an overbarrier resonance [19]. One the other hand, the 
2s -, Id - states are subbarrier Gamov resonances, which can be considered, 
for example, by R - matrix theory [20], or in the framework of the potential 
model with taking into account the continuum [21] into account. Regarding 
to the overbarrier resonances, such resonances do not exist in the complete 
theory as far as we know. At present the quasi-classical approach [19] and the 
complex scaling method [23] are developed for the calculation of characteristics 
of such resonances. Therefore, the determination of the wave function, the 
energy and the width of the overbarrier resonance in the real nuclei with the 
realistic average field is a nontrivial problem. This paper is devoted to find a 
solution of this problem. 



For determination of the resonance energy E^es and the width F of the 
resonant 2p - state, it is necessary to solve the time-dependent Schrodinger 
equation [24]: 



on 



II. FORMALISM 



ih 



dt 



2yU 



A + U{r) ^f(f,t). 



(1) 
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And if we use the substitution ^(r, t) = (f{r)expi^—ijtj, then this partial 
differential equation can be converted into an ordinary one: 



E + ^A- U{r) 



<f(r) = 0. 



(2) 



Taking into account spherical symmetry of the potential U (r) and substi- 
tuting the expression (pnij{r) = Xnij{i")/i" for radial part of wave function, we 
find: 



h'^ (f l{l + 1) 



+ 



2/1 dr^ 2(1 H 



+ U{r) - E, 



nlj 



Xnljir) = 0, 



(3) 



where fi is the nucleon reduced mass, / and j are orbital and total angular 
momenta of the nucleon; n is the principal quantum number, which denotes a 
number of knots in the radial wave function and is connected with a spectro- 
scopic principal quantum number given by = n + 1. For the potential 
U{r) = V{r) + Vs,o. + ycoui, we have the following expressions [16]: 
the Woods - Saxon potential 



V{r) = -Vq h + exp 



Ri 



-1 



a 



(4) 



and the spin - orbital interaction 



where 



Vs.o. = 



r dr 



(5) 



(6) 



The nucleus radius is determined as Rq = tqAs. Here tq is the nucleon radius, 
a is the diffuseness of potential and k = 0.263 fm^ is the parameter for the spin 
- orbital interaction. The depth of potential Vq is calculated in accordance with 
the experimental value of the single-particle binding energy Eiji^d = 4.947Mey 
for the neutron on Ip - shell. 
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Although the Coulomb interaction does not need to be involved within the 
present framework of model description of ^^C, it should participate to the case 
of the proton playing a role. It is given in the form 



1, 



r<Ro 
r>Ro. 



(7) 



The solution of equation (3) will be searched with 

Xnijir) = Unijir) + ivnijir). 



(8) 



As r ^ oo, this solution must satisfy the following boundary condition: 

Xnij{r) — > [--)exp{i{kr - r]ln2kr - — + 6i)}, (9) 



with 6i = argV{l + 1 + 277), where r] is Coulomb parameter. 
From eq (9), we have 



lim = = 

dr h 



and this condition selects the discrete complex value 

Efilj = iEres)nlj 



2^ j- 



(10) 



(11) 



Substituting this expressions for the energy (11) and wave function (8) into 
the radial part of the Schrodinger equation (3), and separating the real and 
imaginary parts, we obtain a set of two inhomogeneuos differential equations: 

+ + ^S.O. + Vcoul. + Vir) - {EresUjl Unijir) = IVnlj ■ Vnljiv), 

-^£^ + + ^S.O. + Vcoul. + Vir) - iEres)nlj\ Vnljir) = -^Klj ■ Unijir). 

(12) 

For Tnij = 0, we have Vnijir) = and Xnijir) = Unijir) = ul^ijir). Then, the 
second equation vanishes and the first one converts into homogeneous type as 



nHii + i) , 



u 



nlj 



r) = Er,ij-v>:i^ir). (13) 
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This homogeneous equation with the Wood-Saxon potential can not be 
solved by the ordinary method. Therefore, it should be solved via the ex- 
pansion on the Sturm - Liouville functions [25] (the detailed description of this 
method is given in Appendix). 

The solution of inhomogeneous differential equation of the second 
order can be expressed through that of the appropriate homogeneous 
equation w^°;^-(r) [25]: 



Unijir) = Cui • u^nijir) + Cu2 ' u^^ijir) f ^} df] 



''nlj 



Here all integrals are indefinte ones. If Fnij = 0, it follows that Unij{r) = u^°ij{r) 
and Vnij{r) = 0. Then, from the first equation (14.1) one obtains, Cui = 1 and 
Cu2 = 0. Similarly, Cyi = Cv2 = follows from the second (14.2). Therefore, 
we have 



Unijir) = wg-(r) 

f if ^nijiO^^:i]im) dm (15.1) 

Vnijir) = -Gu^°i]ir) f ( T UnijiOunijiOd^) dr], (15.2) 

where G = j^l^nlj- 

Substitution of the second equation (15.2) into the first one (15.1) leads to 
the integral equation for Unij (r) . The solution of this integral equation is used 
to determine Vnij{r) by the second equation (15.2). After all, one find 



X 



(16.1) 



Vnij{r) = -Gu^°i]ir 



{fu^:i]iOunij{Od^)dr]. (16.2) 



These integral equations can be solved by means of an iteration procedure 
(of consistent approximations) : the initial approximation is chosen as 

ul°lj(r) — the solution of the Sturm — Liouville problem; 



and, then, the first iteration gives 



4/kr) = 4/ko+A 



where 



X f[u^rdj(T)?dTdvdS,dn; 



At-SW = -Gu^At 



nlj 



(°) 



nlj^ 



The second iteration yields 



v%{r) = v%{r) + A«2].(r); 



where 



AnSJ-(r) 



A4;](r) 



X / ^niji.'^) ■ Au^°i^{r)drdiyd^dr]; 



■ Au^^^iOd^dr^; (21b) 



(17) 



(18) 



(19a) 
(19b) 



(20) 



(21a) 



and etc. 



6 



The boundary conditions (9) and (10) can be rewritten as 



Unij [r) — > - exp [kir + ij— ) sm [hRr - 77 In — - — r - y + O/) , (9a 
Vnij{r) ™ exp {kir + 77^) cos (hnr - r]\n r - y + (5/), (9b) 



I'^iWnijirWijir) + <;j(r)i;„;j(r)]/[w2^/r) +i;^;^.(r)]} = A;/, 
}}^{[unij{r)v'^ij{r) - <z/r)i;nZj(r)]/[u^^/r)+i;^;/r)]} = A;ij. 



(10a) 
(10b) 



In these expressions, k = Ur — ikj and 



kR = ^ cos-, kj = ^sm-, (22) 



2(-E'res)nZj 



P = ^iEres)ll, + r2^^./4, t^V^ = . (23) 



The inverse transformation gives 



{Eres)nlj ^ ^^^e - kj), Tnlj = —ku ■ kj. (24) 

Thus, by having u^°ij{r) and {Eres)nij from the solution of the Sturm - Li- 
ouville problem, the wave functions Uni]{r),Vnij{r) and the width Tnij can be 
calculated by means of the iteration procedure. And, this iteration process must 
converge under the execution of the boundary conditions (9a, b) and (10a,b). 
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III. RESULTS AND DISCUSSION 



The above formahsm has been used for calculation of the wave functions 
Unij{r), Vnij{r)^ the resonance energy {Eres)nij and the width Tnij for 2j9j=i/2 - 
state of neutron in ^^C nucleus. For this case, the Coulomb interaction does 
not play any role, and, thereby, it should disappear from the formalism given 
in the previous section. 

For the Woods - Saxon potential the following parameters were used: 
ro = 1.46 /m, a = 0.67 /m, Vq = 32.3 MeV {Ebind{lpi/2) = -4.947 MeV). 
The wave function u'^ijir) and the energy {Eres)nij were obtained by means of 
solving the Sturm - Liouville problem [16,17] (see Appendix). The number of 
the expansion M in sum (A. 12) was taken equal to 30 so as to reach the enough 
good convergence on energy \E2^^j^ — E2^^~^^\ < 0.1 MeV. The values ^2pf/2 
for the different values M are shown in Fig. 1. The value £"2^1/2 ~ ^-^ MeV 
was obtained at M = 30. The value of the expanding coefficients ^^e 

neutron wave function 1^2^1/2 (^) presented in Table 1. 

For Tnjj = ^2pi/2 = 22.244 MeV, the wave functions W2pi/2(^) and '?^2pi/2(^) 
after 5th iteration were tailed smoothly together with their asymptotical ex- 
pressions, i.e. the boundary conditions (9a,b) and (10a,b) were fulfilled. 

For the time dependent part of wave function, we have 



^f(t) = exp(^-^Enij ■ = exp(^-^{E)nij • • ^^p(^--^ ' ^) = 

= exp(^-i-^(k% - kj) ■ ?j ■ ^^p{^^j^R '^i-t^, (25) 

where the first factor is the wave function oscillating over time, 
and the second factor is the descending function over time. 
There were used relations (11) and (24). 
Notice that the full radial wave function is 

^nij{r.t) = Xnij{r)-^{t). (26) 

If this wave function is considered in the 3 - dimensional space (z, r, t), 
where z = Re^nij (^5 1) oi Im^nij (^5 1) , then one obtains a complicated surface 
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of the second order. 

In the plane {z, r), at t = 0, the Unijir) and Vnij{r) are oscillating and increasing 
in amplitude (i.e. diverging) when r — ^ oo. In the plane at r = 0, the 

real and imaginary parts of Xnijif = 0) • "^{t) are oscillating and decreasing 
in amplitude (i.e. convergent) when t — ^ oo. And if the full wave function is 
considered in the plane (z, t(r)) for t = r/v, where v = is nucleon speed 

at the resonance energy {Eres)nij of 2pi/2 - state [k^ = V^^^^^-^^-)^ _ ^j^g wave 
vector for nucleon in 2pi/2 - state), then 



~, , , , / iki — k]) \ I kn-ki \ , 
^{r,t{r)) = Xnijir) ■ exp\-i — r I • exp\ • r I. (27) 

Such a wave function at r ^ oo is convergent and can be normalized to 1. 
The Fig. 2 shows the real and imaginary parts of the normalized radial wave 
function for the neutron in (2pi/2) - state of ^'^C. With such wave functions, 
the root-mean-square radius (r.m.s.), < r^p^^^ >^ = 11.92 /m, for 2pi/2 - 
state of neutron was calculated. For the neutron in lpi/2 - state, we have 
< rfp^^^ >2 = 3.40 fm. The value of the r.m.s. radius of neutron in 2pi/2 - 
state is much larger than one in lpi/2 -state. Therefore, the neutron 2pi/2 - 
state may be regarded as the "halo - hke" state. 

Generally, the resonance wave function obtained from the time-independent 
Schrodinger equation is exponentially divergent at an asymptotic distance. 
Therefore, it is not easy to calculate the r.m.s. radius with such an unnormaliz- 
able resonance state wave function. One way to calculate the r.m.s. radius is to 
use the complex scaling method [22] or to introduce a convergence factor such 
as exp~"^ provided a ^ oo at the end[23]. It should be noticed that the r.m.s. 
radius calculated by the last method[23], i.e. convergence factor method, is a 
complex number for the resonance state, e.g. < r'^si/2 9.41 + i4:.52fm for 
the overbarrier resonance state 351/2 with Er = 1.63MeV and F = 0.24:6MeV 
(the barrier top is 1.50 MeV). Consequently, the quadrupole moments of nuclei 
will be complex values, too, in their approach. This sequence seems to require 
further consideration and deeper investigation. On the other hand, our method 
of calculation for the r.m.s. radius takes into account the time-dependent part 
of wave function by means of t = fir/hkN. As is explained above, this time- 
dependent part makes the full wave function convergent. Thus, the r.m.s. radius 
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can be calculated without any difficulty. 
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IV. CONCLUSION 



In this paper, we suggested the method of solving the system of inhomo- 
geneous differential equations of the second order with help of preliminary de- 
termined solution of corresponding homogeneous equation. Such an approach 
allows to obtain the wave functions, energy and width for the overbarrier reso- 
nance state. 

This method was used for analyzing 2pi/2 - resonance state of the neutron 
in ^^C nucleus. The value of root-mean-square radius of neutron in this state 
is larger than that in the lpi/2 - state. 

In the framework of this method, we propose to calculate neutron and pro- 
ton resonance 2p - states in ^^C and ^^N nuclei and then to explore them for 
analyzing the electron scattering data and /3 - transition. It can also be a good 
application of our method to analyze the elastic proton scattering from ^^C at 
1 GeV [27]. 

Besides, we suggest in a new fashion to look at the structure of the halo - 
nuclei, in which an existence of the sub - and over - barrier resonant single- 
particle states, as well as the correlated multi-particle states can be displayed. 
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APPENDIX A 

The common theory of the expansion on the Sturm - Liouville functions was 
developed long ago [28,29,30]. It is actually the generalization of the expansion 
in the Fourier series. Examples of application of this theory to the atomic 
and nuclear problems are given in [25]. The Sturm - Liouville equation for 
eigenvalues anij and eigenfunction (fnijir) is 
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^ I ^ - I + ys.o. + Vcoui. + anij ■ V{r) 



(fnijir) = Eoij(fnij{r),{Al) 



where Eqij is the fixed negative number (in our case this is the binding energy 
of nucleon in Ip - shell) . 

The functions iPnij{f) must satisfy the following boundary condition, 

fnijir) ^0 at r ^ and r oo. 



If aoij = 1 and Eqij is the eigenvalue of the Schrodinger equation, then the 
solution of this equation will be the first eigenfunction of the Sturm - Liouville 
problem, i.e. (poij{r) is the ordinary physical wave function. For fixed values / 
and j, the depth of potential will grow with increasing anij ■ If we find such anij, 
that energy of n - state coincides with Eqij, then the obtained wave function 
will be the eigenfunction of the Sturm - Liouville problem. The eigenvalues anij 
are infinite succession of positive discrete numbers, which satisfy the condition: 



and anij — ^ oo at n — ^ oo. It follows from the Hermite conjugate of the Hamil- 
tonian operator for the Sturm - Liouville problem that the functions fnijir) are 
orthogonal with the weight V{r) at the fixed j and E^ij^ i.e. 

I ^n'ijyir)^nijir)dr = -6^^^. (A2) 

Other mathematical aspects of this method (e.g. the problems of sign of 

the eigenvalues, of the completeness and the convergence of the expansion at 

r oo) are considered in works [25] in detail. 

It is very convenient to present the eigenfunctions i-Pnij{f^) in an analytical 
form with parameters calculated by means of modifying quasi-classical method 
[31]: 
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^nij{r) = Nr,ijiSyhxp(^-^^HniS), (A3) 

where N^ij is coefficient of norm, Hn{S) is Hermite polynomial, 5" = dS/dr, 
and S{r) is the correct function to be found by means of numerical calculations. 

As the comparison equation, the equation of the harmonic oscillator type is 
chosen, 

^ + (2n+l-52)$ = 0, (A4) 
and the integral relation is obtained for the definition S{r) [25]: 

where 



pir) = ^ -^{Enij - anijVir)) - Vs.o. - Vcoui. - (^6) 

is the momentum analogous to the quasi-classical one, ri is the least root of 
equation p{r) = 0. 

In our work, we express S{r) in term of 

%) = ^+E^T,(?^^), (.17) 
2 k=i ^ h-a J 



Ak = 



- 5^ 5(rj)cos (^-^ — (^8) 

1 ?=n '^kmax + 2 >^ 



a + b^b — a /(2j + l)7r\ ^^^^ 

2 2 '^^max ~l~ 2^ 

where 7fc(^) is Chebyshev polynomial of k power; a = ri{Enij) and b = 
i"2{Enij) are the turning points, which are defined by the equation: 



Enlj - anlj ■ V{r) - Vs.o. - Vcaul. " ^^^^^^ = 0' (^^0) 
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kmax is maximum number of expansion (in our case, kmax = 30 was took for 
reaching the necessary uniformity) . 

Having such a Sturm - Liouville functions, we can solve the homogeneous 
Schrodinger equation 



+ — o ^ + Vs.o. + Vcoui. + Vir) 



^^:i]ir) = Enij ■ u^:i]ir) (All) 



with the boundary condition Xnij ^ for r — ^ and r — ^ oo, by means of 
expansion u^nij(r) with respect to ^nij{f) '■ 



Kijir) = E dn'ij-^n'ijij')- 



{A12) 



n'=0 



If we substitute this expansion into equation (All), add and subtract anijV{r)^ 
multiply on left-side to Lpnij{T)y[f)^ and carry out integration over r, then the 
following equation can be obtained. 



oo 



n =0 



which is an infinite set of equations. In order to solve this set of equations, it 
is necessary to cut off the summation at some fixed number M. Then, from 
the condition for solving this set of equations, one finds the approximate eigen- 
value eI^j^ = {Eres)nij and the coefficients <i^j, accordingly the eigenfunction 



AtAf^oo, 



rpexact — { fp \ 
^nlj — y^resjnlj- 
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TABLES 



TABLE 1. The values of the expanding coefficients (inii/2 foi" the neutron wave 

function "^2^ 



n 


^nli/2 


n 


dnli/2 


n 


dnli/2 


1 


.67484 


11 


-.02874 


21 


-.00491 


2 


.62089 


12 


.02375 


22 


.00408 


3 


-.29481 


13 


-.01974 


23 


-.00337 


4 


.18248 


14 


.01651 


24 


.00275 


5 


-.12831 


15 


-.01386 


25 


-.00221 


6 


.09306 


16 


.01166 


26 


.00174 


7 


-.06997 


17 


-.00984 


27 


-.00132 


8 


.05489 


18 


.00829 


28 


.00096 


9 


-.04343 


19 


-.00699 


29 


-.00064 


10 


.03521 


20 


.00587 


30 


.00036 
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FIGURES 



FIG. 1. The energy values £^2^1/2 versus the number of expansion, M,for the 
neutron wave function, i.e. M is the number of terms taken in the summation 
(A12). 

FIG. 2. Radial wave function for neutron in 2]9i/2 - state: 
a) real and b) imaginary parts. 
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